The y-sky: diffuse spectral distortions of the cosmic microwave background 



O 

(N 

< 

m 
(N 



o 

u 

43 

Or 

6 

(N 
> 

in 
in 

m 

On 
O 



13 



Cyril Pitrovfl 

Institute of Cosmology and Gravitation, Dennis Sciama Building, 
Burnaby Road, Portsmouth, P01 3FX (United Kingdom) 

Francis BernardearQ 
CEA, IPhT, 91191 Gif-sur-Yvette cedex, (France) 
CNRS, URA-2306, 91191 Gif-sur-Yvette cedex (France) 

Jean-Philippe Uzar@ 

Institut d'Astrophysique de Paris, UMR-7095 du CNRS, 
Universite Paris-VI Pierre et Marie Curie, 98bis bd Arago, F-75014 Paris (France), 
Department of Mathematics and Applied Mathematics, 
University of Cape Town, Rondebosch 7701, Cape Town (South Africa) 

The non-linear evolution of the energy density of the radiation induces spectral distortions of 
the cosmic microwave background both at recombination and during the reionization era. This 
distortion has the same spectral signature as the one produced by the re-scattering of photons by 
non-relativistic hot electrons, the thermal Sunyaev-Zeldovich effect, whose amplitude is quantified 
by a Compton y parameter. A diffuse y-sky is then expected to emerge from mode couplings in 
the non-linear evolution of the cosmological perturbations and to superimpose to the point source 
contributions of galaxy clusters. The equations describing the evolution of the y field and a hierarchy 
governing its angular multipoles are derived from the second order Boltzmann equation. These 
equations are then integrated numerically to obtain the first predicted power spectrum of the diffuse 
y-sky. ft is found to be a remarkable tracer of the reionization history of the Universe. 
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I. INTRODUCTION 



The spectrum of the CMB temperature has been ob- 
served to be extremely close to a Planck [black-body 
(BB)] spectrum According to the standard lore, the 
radiative transfer of the cosmic microwave background 
(CMB) does not induce any spectral distortions so that 
the brightness of the radiation propagating in given di- 
rection n % can be described by a BB spectrum. It follows 
that fluctuations in the brightness map into temperature 
fluctuations. Spectral distortions are mostly associated 
with secondary anisotropies and due to the interaction 
of CMB photons with the electrons of the cosmic plasma 
with a higher temperature or due to electron scattering 
on charged particle. They are respectively referred to 
as thermal Sunyaev-Zel'dovich (SZ) and free-free and arc 
both late time effects 0-0] ■ Any energy injection before 
recombination and at redshifts smaller than ~ 10 6 would 
also reflect itself through a spectral distortion described 
by an effective chemical potential. 

The previous statement is however stricly correct up 
to first order in cosmological perturbations. In general, 
fluctuations in the radiation distribution are induced by 
gravitational effects or through Compton scattering out 



of free charged particles (mainly free electrons). While 
gravitational effects do not introduce spectral distortions 
at any order in perturbations - since the geodesies of pho- 
tons with proportional 4-momcnta are actually the same 
- Compton scattering does generate distortions beyond 
the first order in perturbation theory. The spectral dis- 
tortion arises mostly at the last scattering surface (LSS) 
and during reionization, from the electron flows irrespec- 
tively of their thermal velocity dispersion. We are thus 
dealing with a non-linear kinetic SZ effect, physically dif- 
ferent from a thermal SZ effect. However, as we shall see, 
it has the same spectral signature. 
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The emergence of a diffuse component in the spectral 
distortion is thus a generic feature of the cosmological 
dynamics beyond the linear order. It has so far been ig- 
nored in second order effects treatments. The goal of this 
article is to show explicitly that this spectral distortion is 
characterized by a Compton y parameter, that is similar 
to the termal SZ effect, and then to derive its evolution 
equation from the general second-order Boltzmann equa- 
tion Q (from which we use the results and notation) . We 
then estimate its amplitude and distribution (through its 
angular power spectrum). As we whall discuss, this of- 
fers potentially a new window on the dynamics of the 
cosmological perturbations and on the reionization era. 
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II. DESCRIBING SPECTRAL DISTORTION 



Up to second order in perturbations, we thus have 



Let us first describe the spectral disortion. In the 
Fokkcr-Planck approximation, that is at lowest order 
in the distortions, the distribution function of radia- 
tion 1(E) around a background BB spectrum can be ex- 
panded as @, H 



1(E) = f T (E) + yE~ 3 D [E 3 Df T (E)] 



(1) 



with D = E-Jj^, and where fr(E) is the BB spectrum at 
temperature T which depends only on E/T. 

This decomposition defines the temperature T = T(l + 
O) which is the temperature of the BB which would have 
the same number density of photons, and we can name 
it the occupation number temperature. However, it is 
important to realize that other definitions of the temper- 
ature are possible. For instance, one can define a tem- 
perature out of the brightness I oc J I(E)E 3 dE. On 
the background spacetime, the value of the brightness is 
directly related to the energy density of radiation since 
I — p-y , and we can define the brightness temperature as 
the temperature of the BB which would have the same 
brightness as the actual distribution. One would thus 
define the bolometric temperature Tb by [l7j 



Tb 
T 



1 
I 



(2) 



Indeed both at the background and at first-order level, 
both definitions reduce to the same quantity since in each 
direction of propagation the spectrum is still the one 
of a BB. However this is not the case as soon as colli- 
sions induce spectral distortions, that is at second-order 
in perturbations and beyond. According to the decom- 
position (p}, these two definitions are related through 



T = 



(1 + 4y)4 



The decomposition ([I} defines the y-typc distortion 
which is the only distortion generated by second order 
cosmological peturbations. Such a distortion modifies 
the energy density of the radiation but not the photon 
number density, contrary to the photons produced during 
recombination which conserve none of these densities Q . 
In general the distribution of radiation depends on the 
coordinates (rj,x l ) in space-time and on the direction of 
propagation n l . The temperature of the BB spectrum is 
in principle different in each direction. We will often omit 
the (r],x l ) dependence when no confusion can arise and 
shall express all BB spectra in terms of the background 
BB spectrum of temperature T(rf). The temperature can 
be related T through the fractional temperature pertur- 
bation defined by T{ji,x\n>) = T(r?) [l + Q{r],x i ,vP)]. 
We can then relate the BB spectra through the Taylor 
expansion 



ME) = £ 



-1 



■[ln(l + e)]"D"f T (E). 



(4) 



1 = h - (© - \& 

+ yE~ 3 D[E 3 Dfrp] 



1 



6 2 D 2 /r 



(5) 



In this expansion, the position, time and direction de- 
pendence of 1(E) is to be found in and y. We de- 
compose the directional dependence of in symmetric 
trace-free (STF) tensors (see Refs. @, [Io[ for details) as 



e(n') = e; 



^ O 0i f 77— , where the Qi e 



©u. 



are STF 



and where, conventionally, the lowest multipole (i.e. cor- 
responding to I = 0) is noted 00. Since polarisation is 
also generated at first order, we develop the distribution 
tensor of linear polarisation in electric and magnetic 
type multipoles 



00 

^•(n)=E[ £ ' 

1=2 



kt-2 



.nil 



i m t (iBjjikt 



TT 



where TT means that we extract the trace-free part 
transverse with respect to n % . In practice, we only need 
the polarisation at first order and, since the magnetic 
multipoles are not generated in the absence of primordial 
vector and tensor modes, we can safely neglect them. 



III. EVOLUTION OF THE SPECTRAL 
DISTORTION 

The evolution of y up to second order in perturbations 
can then be deduced directly from the Boltzmann equa- 
tion |2j in which we insert the expansion ©. Includ mg 
polarisation (see Ref. 0, dropping the hatted notation 
to alleviate the equations below), it reduces to 



(3) £^ = 



dy 



-V + M + -^Vijrjn 3 + (0 - 1^(0 - ) 

1 3 1 
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3 

—®ijV k n l n 3 n k - —Qijk'i 



1 i 

3 



11 1) 

—ViVjn^ n JI 



(G) 



-^(vkn* - ®)nW + l-E. ljk nV(v k - fc ) 

where y = y + 2 /2; (ii . . .i n ) denotes the symmetric 
trace-free part of a tensor; r' is the interaction rate, 
given by r' = an e CT where n e is the background density 
of free electrons and ct the Thomson scattering cross 
section; v l is the velocity of baryons in the Poisson (or 
longitudinal) gauge @. Note that only y is a second or- 
der quantity so that all other quantities, appearing in 
quadratic terms, have to be taken at first order only. We 
thus omit the orders of perturbations as no confusion 
can arise. We emphasize that in the r.h.s. r' multiplies 
a second order quantity so that it has to be evaluated 
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at the background level, and the fluctuations of n e , re- 
sponsible e.g. or the Ostriker-Vishniac effect [3], enter 
only at third order in perturbations. Furthermore, y is a 
gauge invariant variable as pointed out in Ref. [7| and as 
can be checked directly from the results of Ref. [Tl[ com- 
bined with the decomposition ([5]). Consequently, both 
the right- and the left-hand sides of Eq. (j6|) are indepen- 
dently gauge invariant. More precisely in this expres- 
sion, the variables 0, 0^, E^ t and v 1 are gauge invari- 
ant quantities which reduce to their usual perturbation 
counterparts in the Poisson gauge. A further simplifica- 
tion arises since y is a pure second order quantity and 
vanishes at the background and first order perturbation 
level: d/drj = d/drj + {dx l / drj)d / dx l + (dn l / drj)d / dn l 
needs only be evaluated at the background level, i.e. as 
d/dr] = d/dr) + n l d/dx l . 

We can extract the moments of the evolution equa- 
tion ([5]) to obtain a hierarchy which would be solvable 
numerically. To simplify the expressions obtained, we 
neglect the multipoles of 6 and E for I > 3. We obtain 



drj 



e+i 

{21 + 3) 



(7) 



where 



c - ^(e i -v i )(e i -v i ) + -e ij e i \ 



(8) 
(9) 



C i 
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va + ™ ( e <* - v a) ( ^> - v i)) + M QikGkj 



10" J 20 



C ijk = -yijk + + \ E {i?j ( Q fe> - v k)) . ( 10 ) 



Cijki — —Vijki + g9(,j9/ci) 



C, 



for £ > 5 



(11) 
(12) 



In the unpolarized case, this result is consistent with 
Eq. (47) of Ref. Q , once the integration over the scatter- 
ing terms is performed. 

In the tight coupling approximation (that is at lowest 
order in powers of k/r') the first order Boltzmann equa- 
tion implies that 0^ oc Qiik/r') 1 ^ 1 for I > 2. In this 
approximation, we can thus discard these quadrupolar 
terms and simplify acccordingly the collision multipoles. 
The tight coupling approximation was shown to be accu- 
rate, at second order in perturbations, during the recom- 
bination fl2l |. During reionization, the same approxima- 
tion can be performed, not because tight coupling holds, 
but mostly because the multipoles 6^, Et and Bt behave 
approximately like jt{krj) whereas v" 1 grows since recom- 
bination and terms quadratic in the baryons velocity thus 
dominate. Again, the spectral distortion considered dur- 
ing reionization is thus a (non-linear) kinetic SZ effect. 



IV. THE DISTORTION HIERARCHY 

For a numerical integration, it proves more convenient 
to use the basis of spherical harmonics rather than the 
STF tensors to develop the angular dependence. When 
working in Fourier space, we choose to align the az- 
imuthal angle with the Fourier mode k. Following the 
convention of Refs. [1, HH, we expand on the basis 



(13) 



with N t = i e y/(2£+ l)/(47r) and obtain our master equa- 
tion 



dyT(k) 



2£ + 3 ' 21-1 

with the collision moments given by 
C (k) = -yo(k), 

(-1) 



T'Ce(k) 



C?(k) 




K l (ki)VL„(k 2 ), 



(14) 
(15) 



n=-l 



C?(k) = -y?(k), 
where we have defined the convolution by 



y m _ n (k 1 )y n (k 2 )|(i6) 

(17) 



and with 



>-/ 


d 3 kxd 3 k 2 r , , 

(2 ; )3 / 2 2 ^(k 1 +k 2 - 


k). 




= Vm (k) - er(k) , 








) 




l(£±m)(l±m- 


hi) 


~ V 2P 




±1 4^7 


j{l±m){l±m- 


1) 


~ V 2P 





(18) 

(19) 
(20) 

(21) 
(22) 



V. THE ANGULAR SPECTRUM OF 
DISTORTIONS 

The spectral distortions are generated only when the 
collisions become efficient enough, that is when baryons 
and photons start to decouple and have different ve- 
locities. For adiabatic initial conditions, the velocities 
of all fluids are initially equal and suppressed on large 
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scales, so that no significant spectral distortion is ex- 
pected on large scales. We can thus safely compute 
the spectrum in the flat-sky approximation, even at the 
reionization time, given that it is reliable beyond I ~ 10 
(see Ref. [l4| for more details on this approximation and 
its domain of validity). The sources S™, or more pre- 
cisely the transfer function from the primordial potential 
$i, on which we integrate along the line of sight are given 
by the implicit definition K. {^(ki, k 2 )<E>i(ki)$i(k 2 )} = 
g(?7)[C™(fc) + y™{k)} where r is the optical depth and 
satisfies dr/dij = —t' with r = 0, and g{rf) = r'e~ r 
is the visibility function. It receives two contributions: 
one which peaks in the middle of the LSS at a comov- 
ing distance which is by definition tlss, and one which 
grows for z < 11, due to the reionization of the universe. 
For the first contribution, the spectrum in the flat-sky 
approximation reads, 

yy 1 f dk 2r dki r dk 1± . 
C\ =-2~ / 7^Zv5 A;isma kl!k P(A;i)P(fc2) 



LSS 



E 



(2tt) 3 

drg(r) £ Q? (k r /k)Sp(k 1: k 2 )e ifc " 



In this expression, we have used the functions 



— -J 



QT(k r /k) = i 



\l-m)\ 
(£ + m)\ 



PP (k r /k) 



(23) 



(24) 



where P (k) is the usual initial power spectrum of the 
primordial gravitational potential $i and P™ (cos /3) are 
the associated Legendre polynomials, in agreement with 
the basis (fT3)) used for the multipoles expansion. The 
Fourier modes k are split into a part k r along the line 
of sight and a part orthogonal to the line of sight, 
the latter satisfying the flat-sky constraint k± = (I + 
l/2)rLSS [HI- ^ is also understood that k 2 = k — ki, and 
ctkj.k is the angle between ki and k. Finally, the distance 
r of an emitting point is defined through r(rf) = r/o — r/. 

The second contribution arising from the reionization 
era can be evaluated using a Limber approximation which 
is a refinement of the flat-sky approximation valid when 
the sources are slowly varying and contribute in a wide 
range of distances, which is our case here. We obtain 



cr 



(2 



irj J drdki r dk 1± \ki r \P(k 1 )P(k 2 ) 



E 



r 



,(25) 



in which it is understood that k r = 0, that is k 2r = —k\ r 
and k± = (£ + l/2)/r. 

The results obtained are reported on Fig. [TJ assuming 
the cosmological parameters of WMAP-5 (l5| with 
instantaneous and full reionization, but allowing the 
optical thickness to vary by a factor 2 compared to 
its WMAP-5 value. The distortion generated during 
reionization strongly dominates over the one arising 
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FIG. 1: The angular power spectrum for the y distortion 
(Cjl y ) induced at recombination (bottom dashed curves) and 
the signal induced at reionization (upper solid curves). We 
report the contributions obtained for r rc ion = 0.043, 0.087 and 
0.174, drawing them with increasing thickness. 



from the LSS. Assuming that between the LSS and 
reionization the baryons are completely decoupled, 
their velocity perturbation scales as v l ~ Qkrj so that 
the ratio between the signal from the reionization and 
the LSS is expected to be of order T re i on (ry re ion/?yLss) 4 - 
Given that 77 rc ion ~ 50//c Gq and ?7lss ~ 3/fc cq , this 
harsh approximation [2(| leads to the rough estimates 
Cf'(reion)/C'f'(LSS) ~ O(10 4 ) in full agreement with 
Fig. [T] Note also that Cj^reion) does not scale linearly 
with Treion, mainly because increasing the reionization 
means that it starts earlier but the higher-r part of the 
integral contributes less given that the integrand roughly 
scales as P 2 v 2 /r 2 ~ (770 — r) A . 

While Cj T = CJ bTb , the fact that T and T B do not 
coincide at second order in perturbations implies that we 
can define two types of bispectra, Bj B ^f^ B and Bf^T^ _ 
At lowest order in spectral distortions, they are related 
by 



TjTTT _ R T B T B T B _ u yTT 



D 



yTT 



<L- (26) 



This indicates clearly the ambiguity of using one temper- 
ature rather than the other one because of spectral distor- 
tions. Fortunately, as we have just shown for the spec- 
trum, the spectral distortions are mostly induced dur- 
ing the reionization epoch during which the terms of the 
form B^^i are suppressed by the fact that the main 
effect in the temperature is coming from the Dopplcr 
effect. Indeed, the Limber approximation tells us that 
most Fourier modes contributions come from the modes 
satisfying k r <C k when the contributions are extended 
in a wide region in time, whereas the Doppler effect is 
more efficient when k r ~ k, and this brings a geometric 
suppression of the coupling of spectral distortions to the 
Dopplcr effect during the reionized epoch. We thus ex- 
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pect the two bispectra to agree numerically, even though 
they do not mathematically 

VI. CONCLUSION 

In this article, we have shown that, beyond the lin- 
ear order in cosmological perturbations, the fluctuations 
in the CMB cannot be entirely described by a tempera- 
ture fluctuation since there is a generic y-type departure 
from a BB spectrum [2l|. We have derived the evolu- 
tion equation of this spectral distortion at second order in 
perturbations from the second order Boltzmann equation 
and we have emphasized that it leads to ambiguities in 
the way the temperature is defined. As a consequence the 
actual temperature bispectrum resulting from the second 
order coupling depends on which temperature is consid- 
ered (in practice though, we found the differences to be 
very weak). 

The spectral distortion is found to be a non-linear ki- 
netic SZ which has the same spectral signature as a ther- 
mal SZ effect (the linear kinetic SZ effect inducing no 
disortion). We also remind that it is different from a 
modulation of the Doppler effect from density and ion- 



isation variation [16j . The amplitude of this disortion 
has been estimated numerically and we have shown that 
the contribution from the rcionization era dominates the 
one arising from the recombination by almost 4 orders of 
magnitude, a consequence of the late time acceleration 
of the electrons. The properties of the y-sky are then 
expected to depend mostly on the physics of reioniza- 
tion, in a way that can turn to be much more precise 
than the CMB polarization [ijj]. However, before draw- 
ing definitive statements, it is necessary to show that such 
effects can be distinguished from the thermal SZ due to 
hot electrons in galaxy cluster halos. This latter effect 
is expected to produce an average signal larger than the 
signal described here. Although tentative to disentan- 
gle the CMB contribution from the galactic emission and 
the unresolved extragalactic radio sources contributions 
of the sky temperature has been made fl9| . it is clear 
that whether those point sources could in the future be 
identified in a large enough number to separate the two 
contributions is still an open issue. 
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different distances are not the same. 
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functions. 



